Free expansion of a Lieb-Liniger gas: Asymptotic form of the wave functions 
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The asymptotic form of the wave functions describing a freely expanding Lieb-Liniger gas is 
derived by using a Fermi-Bose transformation for time-dependent states, and the stationary phase 
approximation. We find that asymptotically the wave functions approach the Tonks-Girardeau (TG) 
structure as they vanish when any two of the particle coordinates coincide. We point out that the 
properties of these asymptotic states can significantly differ from the properties of a TG gas in a 
ground state of an external potential. The dependence of the asymptotic wave function on the initial 
state is discussed. The analysis encompasses a large class of initial conditions, including the ground 
states of a Lieb-Liniger gas in physically realistic external potentials. It is also demonstrated that 
the interaction energy asymptotically decays as a universal power law with time, Ei nt oc t~ 3 . 
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I. INTRODUCTION 

The physics of one-dimensional (ID) Bose gases in 
many aspects differs from the physics encountered in 
higher dimensional systems. For example, the Lieb- 
Liniger (LL) gas of (5-interacting bosons in one spatial 
dimension becomes less ideal as its density decreases 
and eventually approaches the Tonks-Girardeau (TG) 
limit of a gas of "impenetrable-core" bosons [2| as it be- 
comes sufficiently diluted. The interest in these ID sys- 
tems is greatly stimulated by their experimental realiza- 
tion with atoms confined in tight ID atomic wave guides 
[E EL IE < The special features of effectively ID atomic 
gases [6, 0, EI are reflected by properties of nonequilib- 
rium dynamics in these systems, which have become ac- 
cessible experimentally [5[ . The possibility of finding ex- 



act time-dependent solutions for LL [9|, 1C,EI| and TG 
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of particular theoretical interest, as they have the po- 
tential to provide insight beyond various approximation 
schemes. 

Exact solutions for a homogeneous Bose gas with (re- 
pulsive) point-like interactions of arbitrary strength c, 
and periodic boundary conditions, were presented by 
Lieb and Liniger in 1963 [l|. For attractive interac- 
tions, c < 0, exact LL wave functions were analyzed 
in Ref. (23|. The case of box confinement for c > 
was studied in Ref. [2~j] . In the light of recent ex- 
periments [1, 0, HI exact studies of the LL model are 
even more attractive today [IE EE EE EE El]- Be- 
sides providing insight into the physics of ID Bose 
gases, exact solutions can serve as a benchmark for var- 
ious approximations as well as for numerical approaches 
(see, e.g., Refs. [IE E(|). The calculation of corre- 
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lation functions of a LL gas from the wave functions 
is a difficult task; these functions furnish observables 
like the momentum distribution of particles in the gas, 
and were studied by using various approaches (e.g., see 

Refs. HH, EE EE M EE EE EE El EE EE ETEl). 

Time-dependent phenomena in the context of LL gases 
with finite-strength interactions have been addressed by 
using both anal ytical p , [ToL and numerical methods 
(see, e.g., Refs. [43|, [44|])- Irregular dynamics of a LL gas 
was studied numerically in a mesoscopic system in Ref. 
[43| . In Ref. , it was shown that phase imprinting by 
light pulses conserves the so-called cusp condition for the 
LL wave function imposed by the interactions. 

Exact solutions for ID Bose gases are conveniently con- 
structed by using the Fermi-Bose mapping techniques 
[EH, EE Ell- I n I960 Girardeau discovered that the wave 
function of a spinless noninteracting ID Fermi gas can be 
symmetrized such that it describes an impenetrable-core 
ID Bose gas [2|. This mapping is valid for arbitrary ex- 
ternal potentials Q, for time-dependent problems 
and in the context of statistical mechanics (45|. In fact, 
fermion-boson duality in ID exists for arbitrary interac- 
tion strengths (46l. |47|. Furthermore, a time-dependent 
antisymmetric wave function describing a ID system of 
noninteracting fcrmions can be transformed, by using a 
differential Fermi-Bose mapping operator, to an exact 
time-dependent solution for a LL g outlined by 

Gaudin [9]. This method is applicable in the absence 
of external potentials and other boundary conditions. 
Therefore, it is particularly useful to study free expansion 
of LL gases from an initially localized state. 

Free expansion of interacting Bose gases has recently 
attracted considerable attention. It has been utilized in 
experiments to deduce information on the initial state 
(see, e.g., Ref. [48| and references therein), and can 
be considered as a quantum-quench-type problem which 
provides insight into the relaxation of quantum systems 
(see, e.g., Refs. [IE E3| and references therein). Free 
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expansion of a LL gas has been analyzed in Ref. [14| 
by employing the hydrodynamic formalism ||; it was 
shown that the density of the gas does not follow self- 
similar evolution [l4j]- However, in ID Bose systems, 
most exact many-body solutions are given for the TG 
gas 0, EE EE E El- An important result is that 
the momentum distribution of the freely expanding TG 
gas asymptotically ap pro aches the momentum distribu- 
tion of free fermions [H, [l7[- Recently, we have con- 
structed a particular family of exact solutions describing 
a LL gas freely expanding from a localized initial density 
distribution It was shown that for any interaction 
strength, the wave functions asymptotically (as t — ► oo) 
assume TG form. Even though it is generally accepted 
that ID Bose gases become less ideal with decreasing 
density, this intuition is mainly based on the studies of 
a LL gas in equilibrium ground states Thus, a more 
rigorous analysis of the expanding LL gas, which leads to 
more dilute system, but out of equilibrium, is desirable. 
In particular, it is interesting to study the dependence of 
the asymptotic wave functions on the initial state, and 
to see how are the initial conditions imprinted in the 
asymptotic states. 



Here we study the asymptotic form of the wave func- 
tion describing a freely expanding Lieb-Liniger gas, which 
can be constructed via the Fermi-Bose transformation 
and the stationary phase approximation. In SectionHTlwc 
describe the LL model and the Fermi-Bose transforma- 
tion. In Section IlIII we demonstrate that the asymptotic 
wave functions have Tonks- Girardeau structure, that is, 
they vanish when any of the two particle coordinates co- 
incide. The dependence of the asymptotic state on the 
initial state is discussed. We illustrate that the proper- 
ties of the asymptotic wave functions can significantly 
differ from the properties of a TG gas in the ground 
state of some external potential. This study general- 
izes and adds upon our previous result from Ref. [Tlj |. 
as the initial conditions studied here encompass ground 
states for generic external potentials and various inter- 
action strengths. From the next-to-leading order term 
in the asymptotic regime, we deduce that the interac- 
tion energy of the LL gas decays as a universal power 
law in time E mt cx t~ 3 . This is illustrated on a par- 
ticular example in Section |IV[ where we provide further 
analysis of the particular family of time-dependent LL 
wave functions studied in Ref. [111 ]. Explicit expressions 
for the asymptotic form of the single-particle density are 
provided in Section [Vj In Section I VII we calculate the 
asymptotic single-particle density for free expansion of a 
LL gas from an infinitely deep box potential. We compare 
our exact calculation with the hydrodynamic approxima- 
tion introduced in Ref. [8j], and employed in Ref. [14j in 
the context of free expansion, obtaining good agreement 
for all values of the interaction strength. 



II. THE LIEB-LINIGER MODEL 

A system of N identical (^-interacting bosons in 
one spatial dimension is described by the many-body 
Schrodinger equation [l[ 



i=l 



2cS(xi 

l<i<j<N 



Xj)tp B - (1) 



Here, i^b{xi, ■ ■ ■ , XN,t) is the time-dependent wave func- 
tion, and c is the strength of the interaction. It is as- 
sumed that the initial wave function is localized, e.g., by 
the system being trapped within some external potential, 
before, at t = 0, the trap is suddenly switched off and the 
gas starts expanding. We are interested in the behavior 
of i\>b when t — > oo. Here the spatial dimension is infi- 
nite Xj £ (—00,00), i.e., we do not impose any boundary 
conditions. 

Due to the Bose symmetry of the wave function, it is 
sufficient to express it in the fundamental sector of the 
configuration space, R\ : x\ < X2 < ■ ■ ■ < xn, where ipB 
obeys 



.dipB 
'~dt 



N 

£ 

i=l 



(2) 



The ^-interactions create a cusp in the wave function 
when two particles touch. This can be expressed as a 
boundary condition at the borders of i?i [l[: 



f 



d 



dx 



_d_ 



4' B = 0. 



(3) 



X j _|_ 1 — Xj 



These boundary conditions can easily be rewritten for 
any permutation sector. In the TG limit, i.e., for 
c — > 00, the cusp condition implies that the wave 
function vanishes when two particles are in contact: 
1pB(xi, ■ ■ ■ ,Xj,X j+ i, . . . ,x N ,t)\ Xj+1=Xj = [2, HI]. 

Exact solutions of the time-dependent Schrodinger 
equation can be obtained by using a Fermi-Bose map- 
ping operator 0, [ill ] acting on fermionic wave functions: 
If iPf{xi, . . . , xn, t) is an antisymmetric (fermionic) wave 
function, which obeys the Schrodinger equation for a non- 
interacting Fermi gas, 



.dipF 



N 

£ 



dx? 



then the wave function 



where 

o c = n 

l<i<j<N 



sgn^ - Xi 



(-- 


'-)} 


\dxj 


dxij _ 



(4) 



(5) 



(6) 



is the differential Fermi-Bose mapping operator, and J\f c 
is a normalization constant, obeys Eq. (JTJ) Q. For the 
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purpose of completeness we outline, in Appendix [Al the 
proof that the wave function ([5]) obeys both the cusp con- 
dition imposed by the interactions and the Schrodinger 
equation 



III. FREE EXPANSION: ASYMPTOTICS 

In this section we study the asymptotic form of time- 
dependent LL wave functions t/'b.c which are obtained by 
the Fermi-Bose transformation ([5]). All information on 
the initial condition ipB,c(xi, ■ • ■ , xn , t = 0) is contained 
in the initial fermionic wave function ippfei, . . . , xn, t — 



0): 



ip B ,c(xi,...,x N ,0) = j\f c O c ^ F (xi, . . . ,x N ,0). (7) 



The initial bosonic wave function, which can be ex- 
pressed in this way, is assumed to describe a LL gas 
in its ground state when trapped in some external po- 
tential V(x), e.g., in a harmonic oscillator potential, or 
some other trapping potential used in experiments. We 
consider the evolution from this initial state after the 
trapping potential has been suddenly turned off, as stud- 
ied in experiments to deduce information on the initial 
state |48| . The time-dependent fermionic wave function 
iPf{x\, . . . ,XN,t), which freely expands from the initial 
condition ipp(xi, . . . , xjv, 0), can be expressed in terms of 
its Fourier transform, 



(8) 



iPf{xi, . . . , xn, t) — J dki---dkN 

x ip F (ki, . . . , fc^e 1 ^™!^-^)*], 
where uj(k) = k 2 , and 

1 



ip F (ku ■ ■ -,k N ) 



dxi ■ ■ ■ dx 



{2-k) n 

x i/)f{%i, ■ ■ -,XN,t = Q)e l ^i=i K J 



N 



(9) 



By using the Fermi-Bose transformation, the time- 
dependent bosonic wave function describing the freely 
expanding LL gas can be expressed as 



(10) 



?pB,c = I dkf- dk N 

x G(k u fc w ) e J £f=i[ fc ^-^)*], 

where the function G{k\, . . . , k^) is defined as 



G(h, . . . ,k N ) =M c ^F{ki 
x II 

l<i<j<N 



,k N ) 
i) + -(kj -h)}. 



(11) 



It should be noted that G(fc 1; . . . , k^) is not the Fourier 
transform of ipB,c because it depends on Xj through the 
sgn.{xj — Xi) terms. 

The asymptotic form of the wave function (flO]) can 



be obtained by evaluating the integral with the station- 
ary phase approximation. The phase (j) = 5"^-— i \kjXj — 

uj(kj)t] is stationary when dip/dkj = 0. Let {kj} denote 
the fc,-values for which 



dkj 



xj - 2kjt = 0, 



that is, kj — Xj/2t. The phase can be rewritten as 

N 
3=1 

The leading term of the integral in Eq. (fTU|) . as well as 
the next-to-leading term, can be evaluated by expanding 
G(ki, . . . , fcjv) = G({k}) in a Taylor series around the 
stationary phase point {k }: 



^ c =e^« fe '»[G({fc'}) fdk 1 ---dk N e- u ^^- k '^ 



N 

E 



dG({k}) 



{/='} 



dk x ■ ■■dk N (k i ~ k' z ) e - lt ^=^- k ^ 2 



1 ^ d 2 G({k}) 



E 



2! ^ dhdkj 



{k'y 



dk x ■ --dkNih - k' t ){k 3 - k' 3 )e- H ^=i (kl - k ''> 2 + 



(12) 



The remaining integrals in the three terms written out in The third term is nonvanishing only for i = j. Thus 
this expansion can be calculated analytically. The second 
term involving the first derivatives of G({k}) vanishes. 
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Eq. (fl~2|) reduces to 



JV 



TV 



d 2 G({k}) 



{fc'} 



(13) 



From Eq. (|13p we obtain in leading order the asymptotic 
wave function 

-0OO oc i _Ar/2 JJ [sgn(xj - xi) + -(k'j - k[) 

\<i<j<N 

x ^ F (k[,...,k' N )e i ^f=^-^ k >y, (14) 

which is written in a more convenient form in terms of 
the variables £j = Xj/t: 



^«r ff / 2 J] [sgn(0 - 6) + f & -6) 
x ^F(£i/2,...,W2)e* E '=^ 



2c 



(15) 



Equation (| 1 5(1 is the main result of this paper. Evidently 
the asymptotic form of the LL wave function -0oo has TG 
form. Namely, the Fourier transform of a fermionic wave 
function ^p (£i/2, ■ • ■ , is antisymmetric, which im- 

plies that t/'oq is zero whenever £j = (i 7^ j). Further- 
more, is symmetric under the exchange of any two 
coordinates and . This clearly shows that a localized 
LL wave function during free expansion asymptotically 
approaches a wave function with the TG structure. How- 
ever, it should be emphasized that the properties of the 
asymptotic state are not necessarily similar to the wave 
function describing TG gas in equilibrium, in the ground 
state of some external potential. The connection between 
the initial and the asymptotic state is illustrated below. 

In the derivation of Eq. (fT5")) we have analyzed LL 
wave functions which are obtained through the Fermi- 
Bose transformation |5]). This class of wave functions is 
quite general and corresponds to numerous situations of 
practical relevance. Let us discuss the case in which the 
initial bosonic wave function ipso = ipB,c(xi, ■ ■ ■ , %n, 0) 
is a ground state of a repulsive LL gas in an experimen- 
tally realistic external potential V(x), e.g., a harmonic 
oscillator potential. The eigenstates of the LL system in 
free space are of the form 



^ {k} = U{{k})6 c det[e 



ikmXj 1 N 

\m,j=l 1 



(16) 



where the set of N real values {k} = {k m \ m= 1, . . . , N} 
uniquely determines the eigenstate; the normalization 
constant is given by 



1 



AA({fc}) ^ 



{2n) N Nl Yl 



see Ref. [33]. In free space, there are no restrictions 
on the numbers k m . If periodic boundary conditions are 



imposed as in Ref. [l[ (i.e., the system is a ring of length 
L), the wave numbers fc,_must obey a set of coupled 
transcendental equations [l], [2f| [26|, [29|, [39[ which depend 
on the strength of the interaction (see, e.g., Ref. [2i 
The LL eigenstates ^{k\ possess the closure property 
and they are complete [5l|. Thus, our initial state ipso 
can be expressed as a superposition of LL eigenstates, 



4, 



BO 



{k} 

6c^N({k})b({k}) det[e lfcm ^ 
{fc} 



N 



(17) 



where the coefficients b({k}) can be obtained by project- 
ing the initial condition ipBO onto the LL eigenstates. 
By comparing Eqs. J?} and (fT"7]) we find that the initial 
fermionic wave function is 



^fo =Af- 1 J2M({k})b({k})det[ 
{fc} 



\m,j — l " 



(18) 



Since we have assumed that V(x) is an experimentally 
realistic smooth function, also ippo is smooth and differ- 
entiable such that the operator O c can be applied. 

The connection between the asymptotic state (|15l) and 
the initial state tpBo is made through the Fourier trans- 
form of the initial fermionic wave function ipp({k}). 
More insight into the connection between the initial state 
and the asymptotic state can be made by expressing 
ipp({k}) through the coefficients b({k}) utilized in the 
expansion (|TT|> . First, let us note that the coefficients 
b({k}) — b(ki, fe, . . . , fcjy) are antisymmetric with re- 
spect to the interchange of any two arguments fc.; and 
kj (i 7^ j). This follows from the fact that the LL eigen- 
states ip{k} possess the same property, see Ref. (33|. By 
using this property of 6({fc}), Eq. (fl~8"|) can be rewritten 



as 



4=1 K Pj x j 



yj F0 = Af- 1 ]T AA({fc})6({fc}) J^RV^ 
{fc} p 

= kt 1 ^2^2N(kpi,kp 2 , . . . , k PN ) 
p {;=} 

xb(k P1 ,kp2,...,k PN )e i ^ kp ^ 

= Af- 1 7V!^AA({fc})6({fc})e iE f=i^^ . (19) 
{fc} 

By comparing Eqs. (| 19|) and ([8]) we obtain 



^ F ({k})=JV- 1 NW({k})H{k})- 



(20) 



Evidently, the Fourier transform of the initial fermionic 
wave function ipp({k}) is directly proportional to the pro- 
jections b({k}) of the initial bosonic wave function onto 
the LL eigenstates. From this relation we can conclude 
that the asymptotic wave function (|15[) has TG structure 
as a consequence of the antisymmetry of the coefficients 
6({fc}), which originates from the antisymmetry of the 
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FIG. 1: (color online) Contour plots illustrating free expan- 
sion of N — 3 bosons from the ground state of a LL gas in 
a box with infinitely high walls (L = 7r). The left column 
depicts the initial ground state \ipBo(L/2, X2, X3)\ 2 , and the 
right column depicts the asymptotic state i/>oo(0, £2, ?3)| 2 , for 
c = 0.2 (a,b), c = 1 (c,d), c = 2 (e,f), and c = 10 (g,h). The 
density of the asymptotic state is zero when two coordinates 
£i and (i 7^ j) coincide. 



LL eigenstates with respect to kj arguments [33 ] . It is 
also worthy to note that Eq. (|10p . and therefore our main 
result, can be obtained without explicit use of the Fermi- 
Bose transformation; by writing the time dependent LL 
states as tp B ,c = J2{k} &({k})^{fc} exp(-« J2j k]t), and 
after employing the antisymmetry of b({k}) [equivalently 
as in Eq. (fT9]l ] one obtains Eq. (fl~0|) . Formulae (fl"5j) and 
(|20[) provide, under general conditions, the asymptotic 
form of the wave functions for the freely expanding LL 
gas, and the connection between these asymptotic states 
and the initial states. 

For the sake of the clarity of the paper, let us illus- 



trate the asymptotic state of the LL gas on a particu- 
lar example. Suppose that initially the LL gas is in the 
ground state, enclosed in an infinitely deep box of length 
L. The ground state ipso for this potential was found 
by employing the superposition of the Bethe ansatz wave 
functions in Ref. [241] . The coefficients b({k}) can be 
relatively easily found for a few particles by employing 
a computer program for algebraic manipulation (Math- 
ematica). In Fig. [1] we illustrate the initial state and 
the asymptotic state for the case of N = 3 particles, and 
for values of c = 0.2, 1, 2, and 10, by showing the contour 
plots of the probabilities \ipBo{L/2, X2, ^3)| 2 (left column) 
and \tpoc(0, £2, ^3 ) 1 2 (right column). Thus, one particle is 
fixed in the center of the system, while the plots illustrate 
the probability of finding the other two particles in space. 
The left column illustrating the initial states shows that 
the system becomes more correlated with increasing in- 
teraction strenght c and it enters the TG regime for suf- 
ficiently large c (e.g., for c = 10 depicted in Fig. Q] (g) 
the ground state of the system is in the TG regime). The 
right column illustrating the asymptotic state shows that 
the wave function is zero whenever two of the coordinates 
coincide. However, it is important to note that the prop- 
erties of the asymptotic wave functions, even though they 
possess the TG structure, can significantly differ from the 
properties of the TG gas in the equilibrium ground state. 
This can be seen by comparing the asymptotic state in 
Fig. [T](b), and the TG ground state shown in Fig. [l](g). 
The asymptotics of Fig. [1] (b) is obtained after free ex- 
pansion from a weakly interacting ground state (c = 0.2); 
from Fig. [1] (b) we observe that when one particle is fixed 
at zero, there is still a relatively large probability of find- 
ing the other two particles to the left and to the right 
of the fixed one. In contrast, for the TG ground state 
shown in Fig. [1] (g), if one particle is fixed in the center 
of the system, the other two are on the opposite sides 
of that one. Furthermore, by comparing the asymptotic 
states in Figs. [l](b), (d), (f), and (h), we see that their 
properties depend on the interaction strength c. It is 
worthy to mention again that free expansion can be uti- 
lized to deduce information on the initial state (see, e.g., 
Refs. [48| and references therein); free ID expansion can 
distinguish between different initial regimes of the LL gas 
B 

Let us now address the case of attractive interactions. 
For c < 0, the cusp condition assumes a form that is iden- 
tical to that for c > (see, e.g., Ref. [2^]). Therefore, 
by acting on some fermionic time-dependent wave func- 
tion obeying Eq. ^ with the Fermi-Bose transformation 
operator O c <o, one obtains an exact solution for the at- 
tractive time-dependent LL gas in the form O c< q^f (see 
Appendix fXJ); our derivation holds for this family of wave 
functions. Experiments where the attractive quasi- ID 
Bose gas is suddenly released from a trapping potential 
were used to study solitons made of attractively inter- 
acting BEC [52j] ■ Exact studies of such a system within 
the framework of the LL model are expected to provide 
deeper insight into nonequilibrium phenomena beyond 
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the Gross-Pitaevskii mean-field regime, where interest- 
ing dynamical effects can occur |53l. l54j|. 

It should be noted that the time scale it takes for the 
LL system to reach the TG regime depends on the initial 
condition. The next-to-leading term of the asymptotic 
wave function is suppressed relative to the leading term 
by a factor 1/t, as obtained by the stationary phase ex- 
pansion in Eq. (| 13|) . From this we can deduce the scaling 
of the interaction energy, defined as 

E int =2c I dxi ■ --cLxn \ipB,c\ 2 ^2 S(xi - Xj), 

l<i<j<N 

(21) 

as t — > oo. Since the interaction strength c is finite, and 
since the asymptotic density |t/'oo(Clj • • • j Cn, t)\ equals 
zero for any pair of arguments being equal, = £j , one 
concludes that asymptotically the leading term of the 
interaction energy vanishes. Since the first correction to 
the leading TG term of the wave function is of order t -1 , 
and since 8{xi—Xj) — t" 1 8{^i—£ l j), the interaction energy 
asymptotically decays to zero as E- mt oc i~ 3 . This power 
law decay of the interaction energy is illustrated in the 
following section. 



IV. EXAMPLE: FERMIONIC WAVE FUNCTION 
EXPANDING FROM A HARMONIC TRAP 

In Ref. [Tl| . we have constructed a particular family 
of time-dependent wave functions describing a freely ex- 
panding LL gas. The wave functions were obtained by 
acting with the Fermi-Bose mapping operator onto a spe- 
cific time-dependent fermionic wave function, 

r N 2 v v - iv 2 t ^ r xa 1 2 } 

oc a v\- t —^r[t)-— j —^[ 5 ^\ } 

j=i 

xb(t)- N2 / 2 J] [Xi-Xi), (22) 

l<i<j<N 

which describes free expansion of noninteracting fermions 
in one spatial dimension. The initial fermionic wave 
function at t = corresponds to a fermionic ground 
state in a harmonic trap V(x) = i/ 2 x 2 /4 (see, e.g., 
Ref. [Hj]). Here, v corresponds to the trapping frequency, 
b{t) = \/T + t 2 v 2 , and r(t) — arctan(2/£) /V . The limiting 
form of the LL wave function, ipB,c(vib{t), ■ ■ ■ , r/jvb(t), t), 
for t — ► co, was shown to have the following form char- 
acteristic for a TG gas: 

^B,c(m&(*)» • • • , wHt), t) oc b{ty N/2 

i\t2 ■ 2-l N 

x expj-*— r(f) 

x II 9(Vj-Vr) + 0(l/t), (23) 

\<i<j<N 




0.1 1 10 

t 



FIG. 2: (color online) Time-evolution of the interaction en- 
ergy £7int(t), expressed in units of the total energy E. The 
three curves correspond to values of c = 1 (solid line), c = 5 
(dashed line), and c = 10 (dotted line). The straight lines 
depict the asymptotic power law behavior of the interaction 
energy, Ei nt (t) oc t~ 3 (see text for details). 



where g(rf) = \r]\+ivr) 2 /2c. Equation (1 15[1 is a generaliza- 
tion of this result given first in Ref. [ll| • Since Eq. (fT5|) 
was obtained with the help of the stationary phase ap- 
proximation, whereas (|23[) is obtained straightforwardly 
from the exact form of the specific LL wave function (see 
Ref. (HI), ^ is worthy to verify that Eq. (fT5]) reproduces 
Eq. (|2"3")l as a special case. In order to do so, we calculate 
the Fourier transform of the initial fermionic wave func- 
tion, i.e., "4>f(xi, • ■ ■ , xn, 0) from Eq. (|22p . Interestingly, 
the Fourier transform has exactly the same functional 
form as the initial condition in x-space: 

^ F cc e~ ££i k V v [] (kj -hi). (24) 

\<i<j<N 

By plugging this form into Eq. (|15p we obtain: 

V>oo oc t- N ' 2 e-^$/M e WW£ 1 & 

x n [i6--6i + ^-eo 2 ]- (25) 

l<i<j<N ° 

After replacing £j = Xj/t with vr\j = vxj/b(t) which 
asymptotically approaches vrjj ~ Xj/t — £j, we obtain 
the functional form identical to Eq. (|23p . This verifies 
the validity of Eq. ([T5"]) in the special case studied in 
Ref. [HI- 

In order to verify the asymptotic power law decay of 
the interaction energy E- lnt obtained in the previous sec- 
tion, let us calculate the time-evolution of Ei nt for the 
specific family of LL wave functions discussed in this sec- 
tion. We calculate integral (f21~j) for N = 3 particles, and 
v = 2. Given these parameters, E- lnt depends on the 
strength of the interaction c and time t. Figure [3 illus- 
trates time-evolution of the interaction energy for three 
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FIG. 3: (color online) The ratio Ei n t/E as a function of the 
interaction strength c, at three values of time, t = (solid 
line), t = 0.5 (dotted line), and t = 1 (dashed line) (see text 
for details). 



values of c; displayed curves depict the ratio Ei nt (t)/E, 
where E denotes the total energy, which is a constant of 
motion. Evidently, after some initial transient period the 
interaction energy starts its asymptotic power law decay 
Ei n t(t) oc i -3 . It should be noted that the contribution 
of the interaction energy to the total energy depends on 
the interaction strength c. This is illustrated in Fig. [3] 
which shows E- mt (t) /E as a function of c at three points 
in time. At t = 0, the contribution of the interaction 
energy to the total energy is non-monotonous with the 
increase of c; it is zero at c = and in the TG limit 
c — * oo, with a specific maximal value in between. The 
form of the curve is preserved for finite values of t, with 
the evident decay of the interaction energy to zero as 
t — > oo. Note that an equivalent non-monotonous behav- 
ior of the interaction energy as a function of c was found 
for the Lieb-Liniger gas in the ground state for c > and 
with periodic boundary conditions (2o| . 



ASYMPTOTIC SINGLE-PARTICLE 
DENSITY 



Given the asymptotic form of the wave function, 
we finally consider the asymptotic form of the single- 
particle density which is of considerable interest for 
experiment. The single-particle density is defined as 
p c {x,t) = N J dx 2 ■ ■ ■ dx N \i) B ,c(x,x 2 , ■ ■ ■ ,x N ,t)\ 2 . For 
studying asymptotics, it is convenient to define the 
asymptotic form in terms of the rescaled coordinates 
£ - x/t: 



Poo{0 =Noot 



N 



d£ 2 . ..d£Ar|V>oo(£,&, ■ ■ ■ , £,n , t)\ 2 ; 

(26) 



here the normalization constant A/"oo is chosen such that 
J d£ Poc(£,) — N, the total number of particles, while the 
factor t N cancels the trivial time-scaling of the asymp- 
totic single-particle density. 

For the specific asymptotic form of the wave function 
(|2~51) we can analytically calculate the asymptotic form of 
the density for a few particles. As an example, for N = 3, 
the normalization constant is 



oo,N=3 



V2ttV 9 (8c 6 + 48c 4 !/ + 90c 2 v 2 + 45^ 3 ) ' 



while the single-particle density has the following struc- 
ture: 



(28) 



x [32c 6 (3*/ 2 -K 4 ) 

+ 16c 4 (33^ 3 - 3^ 2 £ 2 + £K 4 + C 6 ) 

+ 2c 2 (465^ 4 - 60^ 3 ^ 2 + 90is 2 £ 4 + 20vf + f) 

+ 3^(165^ - 60^ 3 C 2 + 30i/ 2 £ 4 + 4vf + £ 8 )]. 

This expression shows that the Gaussian shape of the 
single-particle density is modulated with the iV-hump 
structure characteristic for the single-particle density of 
a TG gas in the ground state of some external potential. 
The corresponding density (|28|) . in terms of r/ = ^jv is 
shown in Fig. 2 of Ref. It should be noted that such 
an asymptotic form of the single-particle density corre- 
sponds to a particular family of time-dependent wave 
functions obtained in Ref. llj. For different initial con- 



ditions one can obtain a different shape of the asymp- 
totic single-particle density as follows from Eqs. (fT5|) and 
([20)) : the asymptotic single-particle density depends on 
tp F ({k}), that is b({k}). 



VI. COMPARISON WITH THE 
HYDRODYNAMIC APPROXIMATION 

Besides providing insight into the physics of interact- 
ing time-dependent many-body systems, our motivation 
to study exact solutions of such systems is to utilize those 
solutions as a benchmark against various approximations. 
Free expansion of a Lieb-Liniger gas has been studied in 
Ref. |14| by employing the formalism introduced in Ref. 
Q , referred to as the hydrodynamic approximation. This 
formalism can be written in a form of a nonlinear evolu- 
tion equation for a single-particle wave function ipnix, t) 
[see Eq. (9) in Ref. 



,dip H (x,t) d 2 ip H 



dt 



dx 2 



+ V(x)^ H + c 2 f 



(29) 



where \tpH(x,t)\ 2 denotes the single-particle density nor- 
malized to J \ipn(x, t)\ 2 dx — N, while the function/ 
which appears in the nonlinear term is defined in Ref. [8j| , 
and also tabulated in Ref. [19] of Ref. [8j. The potential 
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FIG. 4: (color online) The asymptotic form of the SP den- 
sity obtained exactly (black solid line), and with the hydro- 
dynamic approach (red dotted line). The parameters used in 
the calculation are N — 3, L = n, c = 1 (a), c = 2 (b), c = 5 
(c), and c = 10 (d). (see text for details). 



numerically by solving Eq. (129)) with the standard split- 
step Fourier technique; the nonlinear term in Eq. (|29[) . 
that is, the function f(c/\ipii(x)\ 2 ) is calculated by using 
values tabulated in Ref. [19] of Ref. The asymptotic 
dynamics in the hydrodynamic approach occurs after suf- 
ficiently long propagation, when the SP density starts 
exhibiting self-similar propagation (see also [13]). 

The agreement is qualitatively excellent for all values 
of the interaction strength, and quantitatively excellent 
for c < 1. The width of the SP density as a function 
of £ = x/t indicates the velocity of the expansion of the 
cloud. The asymptotic FWHM (full-width at half max- 
imum) expansion velocity is in good agreement for all 
values of c. The hydrodynamic approximation does not 
reproduce small humps in the SP density, characteristic 
in the TG regime after expansion from the ground state; 
this discrepancy is expected to be smaller if we had cal- 
culated expansion from the ground state with large TV, 
where the hydrodynamic approximation is expected to 
work even better. 



is V(x) — during free expansion. The hydrodynamic 
approximation was used to obtain Eq. (|29[) . which is 
written in units corresponding to the Lieb-Liniger model 
of Eq. (JTJ. The nonlinear equation above reduces to 
the standard Gross-Pitaevskii equation for small inter- 
actions, and to the nonlinear equation from Ref. [56j 
for strong interactions [141 ] . The hydrodynamic approx- 
imation overestimates the coherence in the system, and 
therefore it may not be accurate for analyzing observables 
strongly connected to coherence. However, it is reason- 
able to compare the exact asymptotic form of the single- 
particle density after free expansion with the asymptotic 
form obtained from the hydrodynamic approximation. 

Let us follow upon our example from Section [Mi that 
is, let us consider the asymptotic form of the single par- 
ticle density Poo(£) of a LL gas which is initially in the 
ground state of a box with infinitely high walls; the length 
of the box is L = it. The calculation of the exact SP 
density demands performing multi-dimensional integra- 
tion over N — 1 variables which is not a simple task. For 
this reason, the number of particles in our calculation 
of the exact SP density is N = 3. For the initial con- 
dition of the hydrodynamic approach ipff(x,t — 0) we 
could choose ipH(x,t = 0) = y/N/L within the box, and 
zero otherwise. This would be a good initial condition 
in the thermodynamic limit (large N, N/L — const.). 
However, since for our exact calculation we used N = 3, 
we have chosen, in order to be able to compare be- 
tween the two approaches, the hydrodynamic initial field 
ipH(x,t = 0) = y/n cxact , where n oxact is the exact SP 
density of the initial ground state (this can be calcu- 
lated by employing Ref. [24[). Figure d] displays the 
exact asymptotic form of the SP density, and the hydro- 
dynamic asymptotic SP density. The latter is obtained 



Another possible comparison that can be made with 
the hydrodynamic approximation is the following. The 
LL wave function which is utilized as the initial condition 
in Sec. IIVI and Ref. [ll| is obtained by acting with the 
operator O c onto the fermionic ground state ipFo m the 
harmonic trapping potential V(x) — v 2 x 2 /A. This wave 
function can approximate the ground state only when 
the commutator [O c , V(x)] can be neglected [llj. The 
SP density of this state can be compared with the static 
hydrodynamic density obtained in Ref. [8( for the LL gas 
in a harmonic trap. Due to the properties of the operator 
O c [Til ] and the fermionic ground state in the harmonic 
trap ipFo, it is straightforward to verify that the shape of 
the SP density corresponding to the state O c i/jfo scales 
as p(x) — > p(x/s)/s under the transformation v — > v/s 2 , 
c — > c/s, that is, the shape of the SP density does not 
change under this transformation. The same is true for 
the shape of the (ground-state) SP density obtained with 
the hydrodynamic approach, which has been shown [8J] to 
depend on a single parameter n = (^r)a that is invari- 
ant under the transformation v — > v/s 2 , c — > c/s. This 
is fully analogous to the case of a homogeneous LL gas 
where the only governing parameter 7 = c/n is invariant 
under a simultaneous rescaling of the interaction strength 
c and the linear particle density n [l|. The shape of the 
SP density of the state O c tpFQ (calculated for N = 3) 
agrees with the shape obtained in Ref. Q only in the 
Tonks-Girardeau limit (7/ <C 1) where OctpFo is a good 
approximation for the ground state. If we reduce the in- 
teraction strength c by keeping v fixed, thereby increas- 
ing 77, the two SP densities will no longer have a similar 
shape; this stems from a simple fact that O c ipFo is an ex- 
cited state for sufficiently small values of c, because the 
commutator [O c , V(x)] cannot be neglected, whereas the 
hydrodynamic solution approximates the ground state. 



9 



VII. CONCLUSION 

We have derived the asymptotic form of the wave func- 
tion describing a freely expanding Lieb-Liniger gas. It is 
shown to have Tonks- Girardeau structure [see Eq. [fT5])]. 
that is, the wave functions vanish when any two of the 
particle coordinates coincide. We have pointed out that 
the properties of these asymptotic states can significantly 
differ from the properties of a TG gas in a ground state 
of an external potential [see Fig. [T], The dependence 
of the asymptotic state on the initial state is discussed 
[see Eq. ([20)) ]. The analysis was performed for time- 
dependent Lieb-Liniger wave functions which can be ob- 
tained through the Fermi-Bose transformation ([5]). This 
encompasses initial conditions which correspond to the 
ground state of a repulsive Lieb-Liniger gas in physically 
realistic external potentials. Thus, our analysis charac- 
terizes the free expansion from such a ground state, af- 
ter the potential is suddenly switched off. In deriving 
our main result, Eq. (fl5|) . we have used the stationary 
phase approximation. This generalizes and adds upon 
the result from Ref. jll] which was derived for a particu- 
lar family of time-dependent Lieb-Liniger wave functions. 
We have demonstrated that the interaction energy of the 
freely expanding LL gas asymptotically decays accord- 
ing to a power law, E int oc t~ 3 . Furthermore, we have 
calculated the asymptotic single-particle density for free 
expansion of a LL gas from an infinitely deep box poten- 
tial. We have compared our exact calculation with the 
hydrodynamic approximation introduced in Ref. [8| , and 
employed in Ref. [1J] in the context of free expansion, 
obtaining good agreement for all values of the interac- 
tion strength. As a possible future avenue of research, 
we point out that the methodology employed here for the 
analysis of asymptotic wave functions has the potential 
to be exploited further to study the evolution of various 
observables (e.g., the momentum distribution which was 
studied for a TG gas [l6[) and correlations (e.g., see [57| 
and Refs. therein) during free expansion. 
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In this appendix we outline the proof that the wave 
function ^ obeys both the cusp condition imposed by 



the interactions and Eq. ((2|), i.e., that it obeys Eq. iJTJ. 
Without loss of generality we restrict our discussion to 
the fundamental permutation sector Ri . Let us write the 
differential operator as O c = rii<i<j<Ar ^iji where 



Bij - 



d 



dxi 



(Al) 



We first show that the wave function ([5]) obeys the cusp 
condition ([3|) (see Ref. [33[). Consider an auxiliary wave 
function 



ipAvx(xi, . . .,X N ,t) = B j+ ijO c 1p F 

= i, ,.,<•/•. (A2) 

where the primed operator O'j J+1 = O c /Bjj+i omits the 

factor Bjj + i as compared to O c . The auxiliary function 
can be written as 



ipAux 



d 



d 
dxj 



0'i> F . (A3) 



It is straightforward to verify that the operator 
Bj + i,jBjj + iOj j +1 in front of ipp is invariant under 
the exchange of Xj and Xj+\. On the other hand, 
the fermionic wave function ipp is fully antisymmetric 
with respect to the interchange of Xj and Xj+i. Thus, 
■0aux(£i, . . . , Xj , Xj+i , . . . ,Xn, t) is antisymmetric with 
respect to the interchange of Xj and Xj+i, which leads to 



ipAvx(xi,--.,Xj,Xj + i,...,x N ,t)\ Xj+1 = Xi =0. (A4) 

This is fully equivalent to the cusp condition ([3]), 

Bj+ijipB,c\xj + i=xj = 0- Thus, the wave function flSJ) 
obeys constraint by construction. 

Second, from the commutators [d 2 /dxf, O c ] = and 
[id/dt,6 c ] = follows that if ipp obeys Eq. then 
ipB,c obeys Eq. @, which completes the proof. 

If we use the expression 



sgn(xj 



1 



(- 


--)} 


\dxj 


dxi J 



(A5) 



we obtain O c — rii<i<i<Ar-^w as m Eq. ©, which is 
valid inside any sector of the configuration space (see 
0). Note that for c — > oo, one recovers Girardeau's 
Fermi-Bose mapping [2j, where the operator O c =oo = 
rii<i<7<Ar s S n ( x j ~ x i) maps a noninteracting fermionic 
to a bosonic Tonks-Girardeau wave function. 
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